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We study quantum ferrimagnets in one, two, and three dimensions by using a variety of methods 
and approximations. These include: (i) a treatment based on the spin coherent state path-integral 
formulation of quantum ferrimagnets by taking into account the leading order quantum and thermal 
fluctuations (ii) a field-theoretical (non-linear er-model type) formulation of the special case of 
one-dimensional quantum ferrimagnets at zero temperature (iii) an effective description in terms 
of dimers and quantum rotors, and (iv) a quantum renormalization group study of ferrimagnetic 
Heisenberg chains. Some of the formalism discussed here can be used for a unified treatment of both 
ferromagnets and antiferromagnets in the semiclassical limit. We show that the low (high) energy 
effective Hamiltonian of a (Si, S2) Heisenberg ferrimagnet is a ferromagnetic (antiferromagnetic) 
Heisenberg model. We also study the phase diagram of quantum ferrimagnets in the presence of 
an external magnetic field h (h c i < h < h C 2 ) and show that the low- and the high-field phases 
correspond respectively to the classical Neel and the fully polarized ferromagnetic states. We also 
calculate the transition temperature for the Berezinskii-Kosterlitz-Thouless phase transition in the 
special case of two-dimensional quantum ferrimagnets. 

PACS number: 76.50.+g, 75.50.Gg, 75.10.Jm 



I. INTRODUCTION 

Since the appearance of Haldane's conjecture JjJ in 
1983 the study of quantum antiferromagnetic spin chains 
has been the subject of much research activity in con- 
densed matter physics. According to Haldane, integer- 
spin Heisenberg antiferromagnetic chains have a unique 
disordered ground state with a finite excitation gap, while 
half-integer chains are gapless and critical. The origin 
of the difference between half-integer and integer spin 
chains can be traced back to the topological "0" term in 
the effective non-linear cr-model (NLcrM) description of 
antiferromagnetic spin chains and is believed to be due 
to non-perturbative effects Haldane's original predic- 
tions were based on large-spin (large-.?) arguments and, 
although a general rigorous proof is still lacking, several 
theoeretical developments have helped to clarify the situ- 
ation and there is now strong experimental and numerical 
evidence in support of Haldane's claim (see Refs. [§], j|, 
and H for details and further references). 

Recently, antiferromagnetically coupled mixed-spin 
chains with an alternating array of two kinds of spins 
have also attracted interest among researchers. Inte- 
grable models of mixed-spin antiferromagnetic chains 
were constructed by de Vega and Woynarovich S and 
the simplest case of such chains with spins S = 1 and 
1/2 were subsequently studied Bj. Since these inte- 
grable models are exactly solvable, they are very use- 
ful for studying (quantum) statistical mechanical proper- 
ties. However, they involve complicated interactions and 
may be inadequate for capturing the essential features 



of mixed-spin chains with nearest-neighbor interactions. 
A similar situation arises in the study of Heisenberg an- 
tiferromagnets where Bethe ansatz integrable chains 0J 
with integer spin and no gap seem to contradict Hal- 
dane's claim. These models include special generalized 
interactions and behave quite differently from the pure 
Heisenberg chain. Furthermore, from a more practical 
point of view, it is the mixed-spin chains with nearest- 
neighbor interactions which are supposed to describe real 
ferrimagnetic compounds 0. 

Several researchers have recently carried out analytical 
and numerical studies on the simplest quantum ferrimag- 
netic chains to obtain the low-temperature properties of 
these systems @] . As can be seen through these works, al- 
though ferrimagnetic spin chains exhibit both ferromag- 
netic and antiferromagnetic features, they show some pe- 
culiar, and sometimes surprising, features uncharacteris- 
tic of either the ferromagnet or the antiferromagnet — an 
example being the existence of gapless excitations with 
very small correlation length. It is important to under- 
stand these features more clearly and much theoretical, 
numerical, and experimental work lies ahead. 

In this paper we present a systematic and detailed 
study of the effect of both thermal and quantum fluc- 
tuations in one-, two-, and three-dimensional quantum 
ferrimagnets with nearest neighbor interactions and with 
arbitrary (alternate) spins Si and S2 on the two sub- 
lattices of a bi-partite lattice. The outline of our work 
is as follows: in Sec. II we present a formulation based 
on the spin coherent state representation and discuss the 
role of fluctuations in quantum ferrimagnets. Using the 
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formulation described in Sec. II, we continue in Sec. Ill 
with explicit analytical calculations of magnetization in 
quantum ferrimagnets at finite temperature and discuss 
the low- and high-temperature limits. In Sec. IV we 
re-formulate the quantum ferrimagnetic spin chain as a 
1 + 1 dimensional quantum field theory and explore its 
novel features and (dis)similarities (such as the existence 
of a topological 9 term) with antiferromagnetic Heisen- 
berg spin chains. In Sec. V the physical properties of 
quantum ferrimagnets, as described by the continuum 
NLuM formulation, are examined further and the quan- 
tum phases of dimerized states are discussed in some de- 
tail. The effective Hamiltonian based on a description 
in terms of dimers, and a detailed quantum renormaliza- 
tion group study of ferrimagnetic chains is presented in 
Sec. VI. We will explicitly show that in the low-energy 
limit quantum ferrimagnets can be effectively descibed 
as ferromagnets and a cross-over to an antiferromagnetic 
phase takes place as the temperature is increased. Wc 
conclude in Sec. VII with a summary of our work and 
closing remarks and suggestions for future directions. 



onto those of the dual system at small \K\. This duality 
symmetry is similar to that of the two-dimensional (2D) 
Ising model, as first discussed by Kramers and Wannier 
Q , where a one-to-one correspondence between the high 
and low temperature limits of the partition function can 
be established. In the ferrimagnetic described 
by ([!]), one may think of K as (a fictitious) temperature 
and make an analogy with the familiar 2D Ising model. 
A consequence of this self-duality is that the ferrimag- 
net must exhibit identical behavior in the \K\ = oo and 
K = limits. 

On general grounds, the sign of K specifies the long- 
range (classical) behavior of the system, i.e. the ground 
state has either a non-zero staggered, or a fully polarized 
magnetization for positive or negative -fT's respectively. 
When K > 0, fluctuations around the classical solution 
can be taken into account by reparameterizing the unit 
vectors m; as flQl 



1 - i>iipiZ + ip xi x + ipyiy 



i e A, (2a) 



II. SPIN- WAVE THEORY AND THE 
SEMICLASSICAL LIMIT 

In this section we present a formulation based on the 
semiclassical expansion of the spin path intergral to ob- 
tain the phase diagram of quantum ferrimagnets with 
nearest-neighbor interactions in d < 3 dimensions. 

To obtain the Hamiltonian for the ferrimagnet we be- 
gin by dividing the spin system into two (arbitrary) sub- 
lattices A and B with spins S and K S on each sublattice 
respectively, where \K\ = S2/S1 is a number that spec- 
ifies the relative magnitude of the spins S\(— S) and 
52 (= KS) on the A and B sublattices. Note that the 
sign of K specifies the relative direction of the spins on 
sublattices A and B. The classical Hamiltonian can be 
written as 



lassical 



KJS' 



E : 

<i,j> 



(1) 



where J > 0. We have used the spin coherent state repre- 
sentation to write the spin operators on the A and B sub- 
lattices as Si = Srxii and = KSnii respectively, with 
rrij a 3-component classical unit vector || . As can be eas- 
ily seen, the Hamiltonian ([!]) exhibits a duality symmetry 
between the large and small \K\ limits. One may increase 
the magnitude of spins on sublattice B, or equivalently 
reduce the magnitude of spins on sublattice A. In gen- 
eral the Hamiltonian (|]J) is invariant under the transfor- 
mation K — > 1/K and S — > KS. There is thus a duality 
symmetry between the large and small \K\ values and, 
in this sense, the ferrimagnetic spin system is self-dual: 
The properties of the system at large \K\ can be mapped 



m, = -yl - iwkz + ip xl x + ipyiy ieB, (2b) 

where ipi = ip x i + iip v i is a complex scalar field. The 
macroscopic spontaneous magnetization, , along an 
arbitrary axis, e.g. z, can be taken as the order parameter 
for the ferrimagnet and is given by 



Mf = (S zi ) ieA -sgn{K)(S zi ) 



ieB- 



(3) 



The lowest-lying collective modes of the system can be 
taken into account by linearizing the Hamiltonian in 
terms of the fluctuation fields, ipi 



H R s = H. 



cff — -^classical Bq 



KJS 2 



(v^i + + + Wj), ( 4 ) 



ieA j£B 

where E is the ground state energy of the system. The 
dynamical part of the action consists of the quadratic 
(Wess-Zumino) term || 



S 



wz 



ns_ f hl3 

Wo 



<h ( V fadrlpi -KJ2 \ , (5) 

ieB J 



which reduces to the corresponding Wess-Zumino term 
for the familiar quantum (anti)ferromagnet when K = 
±1. It is important to note that when K < 0, one must 
use Eq. (|2^ ) for the unit vector m,; on the sublattice 
B. This results in a minus sign in front of the intra-site 
interaction terms in the Hamiltonian (^J) . 

The contribution of the low-lying fluctuations to the 
effective action can be calculated by Fourier expanding 
the field ip which yields 
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£E (3 n (k))S n (k) 



k n 



a n (k) 



(6) 



( ._{-ihu n + KJSz KJS^(k) \ 
°^ K )-y KJS 1 {k) iKfiuj n + KJSz J " 1 J 

Here 2 is the coordination number, -f(k) = e~ 4 ^' k , and 

£ is a vector that connects the nearest-neighbor spins, and 
uj n = 2irn/hl3 denotes the bosonic Matsubara frequen- 
cies. The dispersion relation for the collective modes is 
then obtained by evaluating the poles of S n (k): 

m ± (k) = ^^-JSz 



±\jSzJ{K -If + 4K[1- ^| 7 (fc)l 2 ] 



(8) 



and then by applying the analytical continuation of the 
Matsubara frequencies ihu>^ — > zLTiuj^. Note that Eq.(||) 
gives the dispersion relations for K > 0. For K < 0, it is 
easy to find the dispersion relation of the spin waves, by 
exploiting ihuj^ — > fiio^. 

This equation explicitly exhibits the self-duality of the 
quantum ferrimagnet that was mentioned earlier. As 
can be readily seen from Eq. (||), the K — ► 1/K and 
S — > ATS 1 limits yield similar dispersion relations, thus 
verifying the aforementioned self-duality of the quantum 
ferrimagnet in A-space. In general, S n {k) has two posi- 
tive and negative poles depending on the sign of K . For 
K = 1 wc find the classical Neel state with gapless spin 
waves, oc ±fc. In this case both negative and pos- 
itive poles are present. 



[II we will show that 



In Sec. 

the negative poles contribute to the zero point quantum 
fluctuations (ZPQF's) at zero temperature and result in 
a staggered magnetization, in agreement with calcula- 
tions using the Holstein-Primakoff transformation. For 
K = —1, in which case all spins are aligned along the 
same direction, the poles are positive and there are no 
zero point quantum fluctuations. The poles in this case 
consist of gapless (Clt ~ k 2 ) and gapped (0,7 ~ A — fc 2 ) 
modes since the lattice contains a double-basis Brillouin 
zone. Treating K as a continuous parameter (see also 



comments in the paragraph following Eq. for tuning 
phase transitions, we see that at K = (and zero temper- 
ature) the ground state of the system changes abruptly 
and the negative poles change sign. Quantum fluctua- 
tions are responsible for this change and a continuous 
phase transition takes place at K = 0. In the latter, 
the system is completely disjointed, and the spins are 
fully uncorrelated. In this case the system is in the dis- 
ordered phase due to the ZPQF's. For K > there 
are two poles: A negative pole, which behaves like k 2 
as k tends to zero, and a positive pole with a gap. For 



small wave vectors the positive pole approaches the gap 
A(= JSz\l — K\) as k 2 approaches zero. In this case the 
ZPQF's still persist. Since phase transitions at non-zero 
temperature T depend on both T and A, for T <§C A the 
gapped mode can be truncated from the Hilbert space 
and the long-wavelength behavior of the system can be 
effectively obtained from the fc 2 -gapless mode. In this 
case the ferrimagnet mimics the ferromagnet in the long 
distance limit, but in the presence of ZPQF's. In the 
other extreme limit when T ^> A, the gap is negligible 
compared to thermal fluctuations and the system behaves 
like an antiferromagnet. We therefore expect a cross-over 
from the pseudo-ferromagnetic state to the antiferromag- 
netic state at T c (oc |1 — K\). To see the physical reason 
for the presence of the gap more clearly, let us transform 
Eq. (^) to a linearized (space-time) equation of motion 

{ -d 2 -JSz(l- K)d T + 2K(JSz) 2 a 2 W 2 } ip = 0, (9) 

where do is the lattice spacing and r is the imaginary 
time. As we can see, the second term in this equation 
is linear in the time derivative and therefore implies a 
time-reversal symmetry breaking. This term is absent 
when K = 1 (or when T ^> A) which corresponds to the 
antiferromagnet. The significance of this time-reversal 
symmetry breaking term lies in the fact that it can be 
effectively considered as arising from a local magnetic 
field that changes the magnitude of the spins on the B 
sublattice. As a result, according to this picture, the gap 
can be interpreted as the Zeeman energy cost for a spin 
flip. 

In the following section we explicitly calculate the mag- 
netization for quantum ferrimagnets to leading order by 
taking the effect of fluctuations into account. These cal- 
culations further support the above ideas. 



III. MAGNETIZATION 

The effect of fluctuations on magnetization can be com- 
puted by using the (finite temperature) Green's functions 

Gf (t) = Mm <T t ^(t)^(0)), if K>0 (10) 



and 



Gf (r) = Urn (TVVi(T)^(O)), if K < (11) 

where +(— ) labels the sublattice A(B) to which the ith 
site belongs, the ensemble average is denoted by (••■), 
and T T is the time-ordering operator. The magnetization 
can be found from (|^) with (S Z i) given by 



i e A 



(S zi ) = { _ (12) 
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By applying the standard techniques, we can use Eqs. (|7|) 
and (||) to evaluate Gf by integrating out fluctuations 
and summing over the Matsubara frequencies, 



G 



KNS 



E 



KMlt 



KJSz 



nvq - hfi^ 



Km: 



KJSz 



Mil 



(13a) 



g: 



kns^ \ nn+ - Mr k 



r n B (^M7+) + l] 



+ _^__ r[ns(W ) + l]}. (13b) 

It is important to note that Gt is the Green's function 
of the sublattice A for any positive or negative K's, but 
G~ is the Green's function of the sublattice B for K > 0. 
One may easily show that the Green's function of the 
sublattice B for K < can be obtained by Eq.(13t) if 



ns(PhQf) + 1 — > nB((3hflf ). It is easy to show that the 
magnetization (along the z-axis) per unit cell at T = 
satisfies the sum rule: M$ = (S Z i)i^A + (S Z i)ieB = 
\Si — Ss\ when K > 0. In this case, both the magnetiza- 
tion and the staggered magnetization are non-vanishing 
simultaneously. The non-zero M$ changes the low energy 
physics of the ferrimagnets abruptly, i.e., they behave like 
the ferromagnets rather than antiferromagnets. we will 
come to this point in Sec. IV. We therefore are led 



to considering Mo as an additional order parameter for 
the ferrimagnets, i.e., a ferrimagnet with K > is fully 
described by two non-zero order parameters M^ eel and 
Mo. Note that for K < there is no difference between 
our definition for the order parameter and the usual 
magnetization Mq. 

Before proceeding further with the calculation of mag- 
netization for the ferrimagnet, and in order to check the 
consistency of the formalism adopted here, let us first ex- 
amine some well-known results. For K = — 1 we have a 
ferromagnet and the magnetization can be obtained from 
Eqs. @, (H), and (|) which give 



M, 



(14) 



The asymptotic limit of Eq. (14) at low temperature 
(fl£ « T « JSz) can be studied by introducing the 
infrared cutoff A -1 (the typical system size is character- 
ized by A). The second term in Eq. ( p~4|) is irrelevant in 
the long-wavelength limit and we obtain: 



1 



M z - 25 cx < 



-0-02 C (|) (^i 



d/2 



d = 2 (15) 



d= 3 



where £(x) is the Riemann zeta function. We can directly 
see that there are no ZPQF's at T = in any dimen- 
sion and that for any low T at d = 1, 2 the correction to 
magnetization diverges due to severe thermal fluctuations 
and there can be no long-range order. This is in accord 
with the Mermin- Wagner theorem , which is applica- 
ble to one- and two-dimensional systems with a continu- 
ous symmetry. However, in d — 3 the correction due to 
fluctuations is finite and there can be long-range order, 
with the order-disorder phase transition taking place at 
T c ~ zJS 5 / 3 . 

The other familiar case corresponds to K = 1 where 
the classical ground state is the Neel state with hQrt = 
iJSzyT — \jk\ 2 /z 2 . The staggered magnetization in 
this case can be found from Eqs. @, @, and (g) 



M. 



Neel 



2S+\ 



x coth 



N ^ 

h 



PJSz 



| 7 (fe)| : 



l-^l7(fc)| 2 



(16) 



and the correction due to fluctuations at low tempera- 
tures is given by: 



M. 



Neel 



2S-l(x < 



-CK2 



T 

JSz 



log A 



d= 1 

d = 2 



-a 3 



0.006 ( T 
d \ 



JSz) 



(17) 



where = iV _1 5Z fc (l— \ j(k)\ 2 / z 2 )^ 1 ^ 2 in d dimensions. 
As we can see, that the ZPQF's are present at T — and, 
similarly to a ferromagnet, long range order is destroyed 
at d — 1,2 at any temperature; a result which is also 
expected from the Mermin- Wagner theorem. 

The results obtained above for the two cases with 
K = 1 and K = — 1 agree with those that have been 
established previously by using the standard spin-wave 
approximation [Q . (We should like to note that although 
the behavior of (at least) one-dimensional (ID) quantum 
antiferromagnets at zero T requires a treatment that goes 
beyond the spin-wave approximation, the finite T results 
discussed above are valid independently of these modifi- 
cations.) 

We shall now treat K as an adjustable parameter to 
study the general case of the quantum ferrimagnet. The 
fluctuation-corrected (staggered) magnetization can be 
found as before and one obtains: 



AI 



K 



(l + \K\)S + 9(K)-^Y,9( K ^) 



x {n B (f3hn+) - sgn(A>s(/3M^)} 
where 9{K) is the Heviside function and 



(18) 
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g(K,k) 



i 

l+K 



V(l-K)H«(l-|7(fc)/ 2 | 2 ) 



K < 



Jf > 



(19) 



It is easy to check that Eq. ( |l8| ) reduces to Eqs. (|l5|) and 
( |l7| ) when K = —1,1 respectively. To simplify let us 
rewrite Eq. (|l8|) as 



M 



K 



(l + \K\)S + 9(K)-a d -Am d z , 



(20) 



where the identity ns{— \x\) = — 1 — ns(|a:|) has been 
used. Here 



ad 



(21) 



K < 

which contributes to the ZPQF's at T = 0, and 
Am " = ^E^> k)[n B ((3hn+) + n B (f3h\^\)]. (22) 

k 

Note that for K < both fijr are positive and = 0. 

The temperature dependence of the ferrimagnetic or- 
der parameter resides in Amf. We can there- 
fore learn about the behavior of the system by studying 
Eq. ( |22] ) for general T and K. Similarly, the temper- 
ature dependence of the (usual) magnetization Mo (for 
ferrimagnets with K > 0) is: 

Mo = \Si - s 2 \ - -J- Yl M/^+) - n B (/3h\n-\)\ . 



(23) 



As can be seen from this equation, there is no ZPQF 
contribution to the M z at T = 0, and that for any low T 
at d = 1,2 the correction to magnetization diverges due 
to severe thermal fluctuations and there can be no long- 
range order (the Mermin- Wagner theorem flj]), similar 
to ©. 

As has already been mentioned in Sec. ||, for K > the 
compitition between the temperature and the energy gap 
in the system results in different long-range behavior and, 
since A oc |1 — K\, we expect the compitition between T 
and |1 — K\ to determine the order-disorder phase transi- 
tions in the quantum ferrimagnet. For instance, the high 
temperature limit is equivalent to the small |1 — K | (i.e. 
K ~ 1) limit. Let us therefore investigate the asymptotic 
behavior of the staggered magnetization in both the low 
and high temperature limits, for K > 0. ^From now on 
until the end of the paper we concentrate our attention 
to the ferrimagnets with K > 0. 

The low temperature limit. This limit occurs in the 
range |Or| -C T <C |1 — K\ and one may neglect the 
1 — \lk\ 2 /z 2 dependence in Eqs. (||) and ( |2l] ) compared 
to (1 — K) 2 in the k — > limit. As a result, a c i — 2(1 + 
K)/\l — K\ and one obtains a residual contribution to 



ZPQF's with no dependence on the dimensionality of the 

f] 

if |) and |Ofe | = 2KJSzk 2 /\l- 



system. In this limit the action contains two poles, fi^ 



JSz(\l- 
K\ as k 



K\+2Kk 2 /\l- 
-* 0, and we find: 



Amf oc 



^^^A d=l 

IK JSz 11 u L 



^-^logA d = 2 



(24) 



and 



Am d z = 0-02C (^j (1 +K)\1- K\ 3 / 2 



( T 



V 2KJSz 



3/2 



(25) 



when d — 3. Note that the contribution of the gapless 
modes is responsible for the infrared divergences 

in ( pi|) and softening of the long-range order of the sys- 
tem in low-dimensions. The behavior of the system is 
different from an antiferromagnet and the temperature 
dependence of its staggered magnetization is the same as 
that of the ferromagnet (see Eq. jl5|)). This corresponds 
to the pseudo- ferromagnet phase in Fig. |l|. 

The high temperature limit. This limit occurs in the 
range T ^S> |1 — K\ or K ~ 1 and, in contrast to the 
previous case, here we may neglect (1 — K) 2 in Eqs. (||) 
and ( p?i"| ) as compared to 1 — |7^.| 2 /z 2 . [In Sec. IV we will 
explicitly show that at d — 1 the behavior of quantum fer- 
rimagnets in the high temperature (high energy) limit is 
governed by antiferromagnetic exchange interactions. Al- 
though the proof of this statement for dimensions d > 1 
is still lacking, we believe the qualitative picture is the 
same at higher dimensions, where a mean-field descrip- 
tion becomes more reliable.] As expected, in this limit 
the system behaves more like an antiferromagnet with 
the staggered magnetization given by Eq. ([l7]) . We label 
this state as the antiferromagnetic phase in Fig. [|. To 
understand the phase boundary in Fig [j] better, let us re- 
call Eq. (|23]). As can be seen from this equation, at zero 
temperature Mq = when K = 1, but as the tempera- 
ture is increased Mq vanishes at K — K c where K c 7^ 1 
and Mq is suppressed by the thermal fluctuations. Since 
in this region (where Mq = 0), the (staggered) antiferro- 
magnetic order parameter, M Neel , is non- vanishing, the 
system is in an antiferromagnetic phase. From Eq. fl23| ) 
one finds that T oc |1 — K C \S is the temperature at which 
a cross-over from the pseudo-ferromagnetic phase into 
the antiferromagnetic phase takes place. 

The order-disorder phase boundaries can now be iden- 
tified in terms of the staggered magnetization M^ . For 
d = 3, and far enough from K = 1, the critical tempera- 
ture is given by 



T c oc JSz\K\l\l-K 



1/3 



(26) 



Note that close to the K = 1 limit the critical tempera- 
ture can be obtained from Eq. ( |l7| ) by solving M^ eel = 0. 
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Fig. [l] shows the qualitative phase diagram for the quan- 
tum ferrimagnets. The phases denoted by antiferromag- 
net and pseudo-ferromagnet are separated by T c ~ |1— K\ 
as discussed above. Due to the Mermin- Wagner theorem, 
and as has been explicitly calculated here, there can be no 
long-range order at non-zero temperatures when d = 1, 2 
and hence no phase transitions either. It is noteworthy to 
cmphesizc that at d = 2 the staggered magnetization and 
the usual magnetization are diverged logarithmically, as 
one can see from Eq.(|2^). This is the reminiscent of the 
low-temperature quasi-long range order phases, where 
the vortices bind in pairs. We therefore expect a non-zero 
temperature Berezinskii-Kosterlitz-Thouless phase tran- 
sition [[l3) takes place when the quasi-long range order of 
the 2-dimensional ferrimagnets is broken (when the vor- 
tices are deconfmed). We will come to this point in Sec. 
0. At K = 1 the ferrimagnetic gap (A cx |1 — K\) col- 
lapses and the Neel state appears. The phases of the sys- 
tem now depend on S, i.e. for half-integer S the model is 
critical and the correlation function of the staggered mag- 
netization has a power law decay, while for integer S the 
spin chain is disordered with a Haldane gap. This com- 
pletes our study of quantum ferrimagnets at zero- and 
non-zero temperature within the semiclassical framework 
of spin- wave approximation. Our calculations include the 
effects of both quantum and thermal fluctuations to lead- 
ing order. 



IV. THE CONTINUUM LIMIT 

In this section we re-formulate the K > and the zero- 
temperature ID quantum ferrimagnetic chain defined by 
(|I|) as a quantum field theory in 1 + 1 dimensions. As 
will be seen below, the continuum action obtained here 
is a generalization of the familiar (NLcrM) formulation 
of the ID Heisenbcrg antifcrromagnct. There are many 
advantages in using continuum field theoretical descrip- 
tions of (quantum) spin systems. These include the pos- 
sibility of studying non-perturbative aspects of the model 
and making the existence (or lack of) certain symmetries 
manifest. A notable example is, of course, the original 
mapping of the ID Heisenberg antiferromagnet onto the 
NLcrM which led Haldane to predict that integer-spin 
chains have a finite spin excitation gap in their spectrum 
and half-integer spin chains do not. This difference be- 
tween integer and half-integer spin chains is understood 
to be due to non-perturbative (topological) effects and 
could not have been obtained by using the standard spin- 
wave theory alone. 

As before, we use the spin coherent state representa- 
tion to study the ferrimagnetic chain in the continuum 
limit. Assuming a total number of 2N sites and periodic 
boundary conditions, it is straightforward to obtain the 
action 



S[m] = S S wz [m,] - KS £ S wz [m,] 
+ dr KJS 2 V] m.i ■ rn, 

JO w„-^ 



<ij> 

iST A - iKST B + H classical j 



(27) 



where ^a(b) is the (topological) Berry phase associated 
with the Neel field on the sublattice A(B) and 

Swz = dt drm(t, t) • [d T m(t, r) x d t m(f, t)], 
Jo Jo 

(28) 

is the Wess-Zumino term. The action S is invariant un- 
der the transformations K — > 1/K, S — ► KS and the 
interchanging the sublattices A and B. (Note that the 
Wess-Zumino term has already been introduced in (|^), 
albeit in a different form. The formalism used in Sec. I, 
in which complex scalar fields ipi were used, is akin to the 
Schwinger boson representation and is very useful at both 
zero and finite temperature. However, for our present 
(T = 0) purposes it is most convenient to use the formal- 
ism adopted here.) The significance of the Wess-Zumino 
term in the continuum treatment of the Heisenberg an- 
tiferromagnet, i.e. when K = 1, can be understood as 
follows ||: On a bipartite lattice, and starting with a 
(classical) Neel state around which the fluctuations are to 
be taken into account, the staggered configuration leads 
to a change of sign of the exchange term in the action ( |27| ) 
to that of a ferromagnet. It is thus the staggered sum of 
the Wess-Zumino terms, a purely quantum mechanical 
contribution, which distinguishes ferromagnets from an- 
tiferromagnets. In order to take the effect of fluctuations 
into account one writes as the sum of a slowly vary- 
ing part, the Neel field n.;, which represents the order pa- 
rameter, and a small rapidly varying part, the transverse 
canting field 1^ , which roughly represents the average spin 

il 



Hi; = ni + (~iy a u, 



(29) 



with the constraints n? = 1 and - 1^ = 0. In the targe- 
s' limit one obtains the well-known Haldane's mapping 
Q of the Heisenberg antiferromagnet onto the (effective) 
0(3) non-linear er-model by integrating out the canting 
fields (the fast modes). The resulting action consists of a 
NLcrM term through the (Euclidean) Lagrangian density 



£nl<tM ( n ) 



2.9 



-(d T n) 2 +v s (d x nf 



(30) 



where g = 2/5, v s = 2a JS, as well as a topological 
surface term expressed in terms of the two-dimensional 
Pontryagin density of the hedgehog fields 



Aop(n) 



8tt 



ey-n- (diii x dju), 



(31) 
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where 9 = 2nS, and £y is the two-dimensional antisym- 
metric Levi-Civita tensor. This term results in a purely 
topological contribution to the action through 6W, where 
W (the winding number) is an integer which appears 
through e lSW in the path integral. (Note that the topo- 
logical term ( |3l| ) is absent when d > 1 fil|.) Since 9 
appears in the path integral through e l9 , it matters 
only modulo 2-7T, and therefore 9 — 0, tt for integer and 
half-integer spin chains respectively. Haldane conjec- 
tured that all half-integer (integer) spin chains are gapless 
(gapped) jlj. It has been known for sometime (T^| that 
the NLctM with 9 = is always disordered and has ex- 
ponentially decaying correlations in the long-wavelength 
(strong coupling) limit. At 9 = it, since the spin- 1/2 
chain is known to be gapless from Bethe's exact solution, 
one expects algebraically decaying correlations, provided 
that the (large- S) mapping onto the sigma model remains 
valid down to S = 1/2. 

We will now examine the more general model with K ^ 
1 which corresponds to the ferrimagnetic chain. In this 
case the action ( p7j) can be rewritten as 

S[n] = -—^—iS(YA + Tg) 
1 + K 

classical? 

(32) 

and in the continuum limit we obtain the Euclidean ac- 
tion: 



C 



r n/3 r 
S[n] = J dr J dxC(n), 



(33) 



where 



£(n) = iM A(n) • <9 T n + £ N La M (n) + £ top (n). (34) 

The first term in Eq. (|34| ) is the usual (dynamical) Berry's 
phase of a quantum ferromagnet, with Mq = |1 — K\S/clq 
the magnetization per unit cell (pair of sites). It is this 
term that results in the ferromagnetic branch of the spin 
waves and corresponds to the trajectory of spin over a 
closed orbit on the unit sphere in the presence of a unit 
magnetic monopole at the center. The contribution of 
the first term in C(n) is equivalent to the area enclosed 
by this trajectory and since either of the smaller or the 
larger enclosed areas on the unit sphere must lead to the 
same Berry's phase, the magnetic moment per unit cell, 
i.e. 2Moao, must be quantized with integral values p7| . 
A corollary which follows from this last observation is 
that the parameter K is also quantized. It is interesting 
to note that although at the spin- wave level one can con- 
sider the parameter if as a continuous parameter, when 
including topologically nontrivial spin histories, the the- 
ory is inconsistent unless K is quantized. [Note, however, 
that in the large-5* limit K can be considered as nearly 
continuous as long as S is large on both sublattices.] The 
remaining terms in Eq. (p4j) are the NLctM 



with 



NLctM = 



g' = 



2.9' 



(K + 1)5' 



-(9 T n) 2 +^(^n) 2 



Aa Q JKS 



and the topological term 



C 



.K+l 



top 



Sn • (<9 r n x d x n). 



(35) 



(36) 



(37) 



As in the uniform ID Heiscnbcrg antifcrromagnet, the 
topological term is proportional to the Pontryagin den- 
sity of the hedgehog fields in two dimensions and, simi- 
larly to the antiferromagnet, it is straightforward to show 
that for d > 1 the topological term is identically zero. 
The equations of motion are readily obtained from (|34| ) 
which gives 



M (n x d T n) a + -L (. 

9 



-dl + n ■ <9 2 n) n Q 



n • V 2 n) n a = 0, (38) 



where a Lagrange multiplier has been used to enforce the 
constraint n(r) ■ n(r) = 1. Using the linearized approxi- 
mation, it is then straightforward to obtain the spin- wave 
dispersion relations 



± iJK\l - K\S 




(1 + Kf 
(1-Kf 



k 2 ± 1 



(39) 



which agrees well with the dispersion relations obtained 
from spin-wave theory. 

An interesting new feature of the NLerM-type formula- 
tion of the ferrimagnetic chain is a topological term with 
9 = (K + l)irS, which gives the familiar result 9 — 2ttS 
when K = 1. As can be easily seen, 9 is always equal to 
one of 0, 7r, or ±7r/2 (mod 2ir) for ferrimagnetic Heisen- 
berg chains. A question that immediately comes to mind 
is whether (in analogy with the Heisenberg antiferromag- 
net) the topological term for the ferrimagnetic chain leads 
to markedly different behavior for different values of 9 or 
not. The answer is that it does not, and this is due 
to the fact that the ultimate long-wavelength behavior 
of the ferrimagnetic chain is governed by the ferromag- 
netic term in ([54]) and the presence of the topological 
term leads to no surprises in this case. This is consis- 
tent with fact that the zero-temperature ground state of 
the Heisenberg ferrimagnet is ordered and does indeed 
correspond to a state with staggered spins and constant 
long-range spin-spin correlations fa]. One may show this 
explicitly by making use of a renormalization group ap- 
proach for action (|33|), i.e., the ferrimagnets have a stable 
ferromagnetic fixed point. However, (periodically) alter- 
nating quantum spin chains other than the Heisenberg 
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ferrimagnet may possess a topological term which can 
dramatically affect the long-wavelength behavior of the 
spin chain due to (strong) quantum fluctuations at zero 
temperature. An example which illustrates this type of 
behavior is the model studied in |l8| which consists of 
a periodic array of two kinds of spins. Models such as 
those studied in ]Is[], where the topological 9 term does 
indeed play a role, are likely to lead to novel interest- 
ing behavior in quantum magnetic systems with useful 
practical applications. 



V. QUANTUM PHASES 

The semiclassical methods that we have used in the 
previous sections to explore the physics of quantum fer- 
rimagnets are, generally speaking, effective descriptions 
which are based on approximations in the continuum 
limit. However, although these powerful techniques have 
enabled us to learn much, to understand the existence 
and cross-over of (quantum) phases in quantum ferri- 
magnets other approximations need to be made. In this 
section we will present a (continuum) formulation based 
on dimerization which is paricularly powerful and well- 
suited for a description of quantum phases. 

The energy eigenstates of each dimer are specified 
through the total angular momentum, I = \S\ — Si\,\S\ — 
S2 1 + 1, • • • , Si + S2, and the ground state and the first 
excited state of a dimer separated by a gap proportional 
to the exchange coupling J. In contrast to a dimer 
with identical spins, whose ground state is unique, the 
ground state of a dimer with alternate spins is degener- 
ate. This degeneracy of the ground state dramatically 
alters the low energy physics of the ferrimagnets; making 
them behave more like ferromagnets rather than anti- 
ferromagnets. However, turning on the interaction be- 
tween the dimers removes the degeneracy (allowing the 
localized states to hop to neighboring dimers) and a gap- 
less band for the lowest I appears. This corresponds to 
the gapless spin wave mode (as obtained in the proceed- 
ing sections) for the intra-band transitions with the dis- 
persion law Wfc ~ k 2 . On the other hand, for isolated 
(non-interacting) dimers the gap between the lowest an- 
gular momentum (£ = \S\ — S2D and the first excited 
state (£ + 1) is robust and persists even for an alternate 
spin chain with a uniform interaction between the dimers. 
This corresponds to the gapped mode of the spin waves 
with inter-band transitions. This gap, which is propor- 
tional to J\Si — S2I, results from transitions between dif- 
ferent bands with different angular momenta £. We will 
return to this point in the next section (Sec. |Vl| ) where 
we show that the ground state of one-dimensional quan- 
tum ferrimagnets is always gapless and ferromagnetic, 
no matter what the values of S\ and £2 are, as long as 
Si + S 2 (K + 1). 



The phases of the system at non-zero temperatures 
can be understood as follows: At low temperature the 
phase of the system can be characterized by gapless tran- 
sitions between the states inside the first band and the 
system behaves like a ferromagnet. At high temperature 
the phase of the system is characterized by gapped tran- 
sitions from the lower to the higher bands {I — > I + 1) 
and the system exhibits antiferromagnetic features in this 
phase (see Fig. [l]) . (Note that due to the Mermin- Wagner 
theorem there can be no long-range order in ferrimag- 
netic spin chains in either of these cases when d = 1,2, 
although at d = 2 the system exhibits quasi-long range 
order.) 

As we have seen before, the ferrimagnetic spin-waves 
consist of both gapless and gapped modes. In the low 
temperature limit the latter mode is irrelevant. At d = 1 
and T = the low (high) energy physics of quantum fer- 
rimagnets is effectively like that of a ferromagnet (anti- 
ferromagnet) which is formed by the chain of (dimerized) 
unit cells with magnetic moment Mq(= \Si — S^D- This 
will be shown explicitly by a renormalization group cal- 
culation in the following section. Applying an external 
magnetic field, h, leads to gapped spin- waves; in which 
case the energy cost for ferromagnetic transitions is pro- 
portional to the Zeeman splitting factor. In fact, the 
dispersion relation for spin waves based on the classical 
ground state of the ferrimagnets (classical Neel state) can 
be obtained from Eq. <M) 



-hS 



\K-1\ 



tJ s * 



+ \jSzJ{K - l)^ + AK[l - 4l7(fc)l 2 ]- (40) 



As one can see, the effect of the external magnetic field 
is to suppress the ferrimagnet gap. Therefore the ground 
state of the ferrimagnet corresponds to the staggered 
configuration of spins on sublattices A and B, unless 
h c i > 2Jz\K — 1|. At this point the staggered state be- 
comes unstable against the canted phase (spin-flop state) 
of the spins when the spectrum becomes soft at k = 0. 
When the external magnetic field exceeds h C 2, the sys- 
tem will be in a saturated ferromagnetic phase, with a 
quantized magnetization per unit cell. This is shown in 
Fig. [|. The h C 2 = 2Jz(l + K) is obtained by using the 
dispersion relation of the spin waves based on the fully 
polarized state of the ferrimagnets. It is the lower-bound 
of the external magnetic field, in the sense that the spin 
waves become soft at k = 0. Clearly, this reveals similar- 
ities between the ferrimagnets and the XX Z + h models 
p2[ . However, we will not pursue this interesting con- 
nection any further and leave it for future investigations. 

As has been shown here, the problem of a d- 
dimensional quantum ferrimagnet can be mapped onto 
a classical d + 1-dimensional non-linear cr-model with or 
without a topological term for spin systems in d = 1 or 
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d 7^ 1 respectively. Turning back to Eq. (p4[) , let us now 
reparameterize the (3-component) unit verctor n as: 

n(r) = (sin#(r) cos y(r), sin#(r) sin<p(r), cos 0(r)) , (41) 

where r denotes the position of n in d + 1 dimensions. 
On the other hand, as we have seen previously through 
the spin wave calculations, when h c \ < h < h C 2 the ferri- 
magnet is in the canted phase, which contains an in-plane 
gapless mode — even in the presence of an external per- 
pendicular magnetic field — and this gapless mode can 
significantly affect the low energy physics JTgj— |2l|] . In 
this connection, it is convenient to think of 6{r) as a 'gap 
field variable' when the 0(3) symmetry is broken. We 
can then integrate out the gapped mode and obtain an 
effective Lagrangian for the in-plane gapless mode (on a 
c?-dimensional lattice) : 



2v 's V- / \ 
-— cos{tpi-tpj). 

<i,j> 



(42) 



The first term in (ji^) is the usual ferromagnetic Berry's 
phase and the last term corresponds to the classical XY 
model. In the continuum limit, the XY model is sim- 
ply the Laplacian of the field ip(r) and by introducing 
momenta conjugate to ipj 



dC 



1 



d{dtfj) g'v' s 



one obtains the effective Hamiltonian 



(43) 



1 d 



i dipj 



2vi 



fJ 



f cositpt -Cpj), 



(44) 



<i,3> 



where we have used the quantum mechanical representa- 
tion of the number of the flipped magnetic moments per 
unit cell 



\Mj 



oc 



difje 



iipj Mj 



(45) 



which is canonically conjugate to the classical in-plane 
orientation [M . The Hamiltonian ([l4|) is the boson Hub- 
bard model (the quantum rotor with a non-zero minimum 
angular momentum). The superfluid (Mott-insulating) 
phase of bosons can be obtained when j' < 2 (g' ^> 2). 
When g' — > oo, and for any value of Mo, except odd 
integr 2Mo's, the ground state consists of the gapped 
Mott-insulating phase, and the classical ordered phase is 
destroyed by quantum fluctuations. It then follows that 
at T — 0, a ferrimagnet with arbitrary (but unequal) 



spins Si and S2 is actually in a superfluid phase of the 
condensated bosons. Within this description we can also 
see that at d — 2 quantum ferrimagnets exhibit critical 
behavior where a Berezinskii-Kosterlitz-Thouless phase 
transition |l3| takes place at Tbkt — nv' s /g'. Further- 
more, the magnetization M (= (Si) + (82)) does not have 
a quantized value and varies continuously as a function of 
the couplings g' and the external magnetic field, h. This 
is the partially magnetized phase and/or the spin-flop 
phase. 



VI. THE QUANTUM RENORMALIZATION 
GROUP APPROACH 

As has been discused in the previous section, it is pos- 
sible to capture some important features of the physics 
of quantum ferrimagnets through effective (field) theories 
based on dimerization within the semiclassical approxi- 
mation. In this section we wish to show that an effective 
theory based on dimerization and not necessarily reliant 
on some underlying assumptions (for example the large- 
S limit) in semiclassical treatments can also reveal much 
about the behavior of these systems. In particular, we 
will show that a quantum renormalization group (QRG) 
study of ferrimagnetic quantum spin chains yields good 
quantitative and qualitative results. As opposed to other 
powerful techniques, such as the density matrix renor- 
malization group (DMRG) p3| , the QRG approach is 
much less complicated and intuitively straightforward. 
The QRG is based on the much familiar standard block 
renormalization and is also analytically quite tractable 
p5[ (see also [^6j for a modern treatment and [M for 
some recent applications). 

To apply the QRG technique one begins by dividing up 
the spin lattice into small blocks (see Fig. ^) and obtains 
the lowest-energy states (\a)) of each isolated block for 
a particular boundary condition — for example, for the 
dimerized states that we will use here open boundary 
conditions for isolated blocks are chosen. The effect of 
inter-block interactions is then taken into account by con- 
structing an effective Hamiltonian H e ^' which now acts 
on a smaller Hilbert space H e ^ embedded in the orig- 
inal one. In this new Hilbert space, each of the former 
blocks is treated as a single site and the technical way 
of implementing this idea is to construct an embedding 
operator 



eff 



H 



and a truncation operator 

Q f : H — -> H eff 
by demanding the commutativity of the following dia- 
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gram pfl): 

H eff H 

H £ ff | i H 

i.e. such that: QH e ^ — HQ. ^From the last relation 
one obtains the effective Hamiltonian as : 



H eff = Q^HQ. 



(46) 



Note that the operators Q and satisfy the relation 
Q^Q = ^-H'ff but QQ^ ^ In- More precisely, one divides 
H as H = H B + H BB (see Fig. |) where H B = J2i hf is 
the sum of block Hamiltonians and H BB — Yl<i j> hfj 
is the sum of inter-block Hamiltonians. Then Q — Y[j Qi 
with Qi given by 



Qi 



E 



a {a 



(47) 



where m is the number of low energy states that are 
kept. Note that since each Qi acts trivially on all the 
other blocks J ^ I, [Qi, Qj] = 0. 

Let us now specialize this construction to a one- 
dimensional quantum ferrimagnet with Si = 1/2 and 
S2 = 1- Here, each cell labeled by i is considered as an 
isolated block in the renormalization group procedure. 
The block Hamiltonian (H B ) and the inter-block Hamil- 
tonan (H BB ) can be written as 



N 



H 1 



>2,i) 



(48) 



(51) 



where a stands for +, — , z and S' is the matrix represen- 
tation of the spin 1/2 operator in the effective Hilbert 
space of the block-renormalized chain. The effective low- 
energy Hamiltonian {H^J) can be obtained from 

K" = Q i H 2 N(J, S)Q = Q^H B Q + Q^H BB Q. (52) 

Using S 2j4 • S M+1 = i(S + 2 ,iS- M+1 + S~ 2 ,iS + i )i+ i) + 
S z 2 jS z i along with Eq.(51) one has 



t t — 4 / / 



from which it follows that 



N 



Qi H BB Q = ^J6^2S' i -S' i+1 . 
J i=i 

^From this last equation and the equation 
Q^H B Q = -NJ, 



(53) 



(54) 



(55) 



we see that the effective low-energy Hamiltonian 
(HlJwi 1 / 2 , 1)) describing a {Si = 1/2, & = 1) ferri- 
magnetic Heisenberg chain is the same as that of a spin 
1/2 Heisenberg ferromagnet with a renormalized coupling 

constant, j/ ^ 2,1) = (4/9) J5 > 0, as given by 



N 



<i(l/2,l) 



H 



BB 



N 



Si 



(49) 



where we also assume that the spin chain obeys the pe- 
riodic boundary conditions: Si( 2 ).at+i = Si( 2 ),i- The 
Hilbert space of each isolated block consists of a 5 = 3/2 
and 5 = 1/2 multiplet corresponding to (J/2) and {—J) 
in energy. We choose the 5 = 1/2 multiplet base kets as 
the kept states in the QRG procedure to find the effec- 
tive low-energy Hamiltonian. (The higher energy states 
of the 5 = 3/2 multiplet are neglected.) If we specify 
these states by |5, M), where 5 is the total spin and M 
is its z-component, we can construct the embedding op- 
erator of each block as 



Qi 



I 1 -,- 1 - 

'2' 2 



(50) 



where |+) and |— ) are the renamed base kets in the ef- 
fective Hilbert space of each block. 

By using Eq.(^) one can find the effective spin oper- 
ators of the renormalized chain: 



This result is in complete agreement with the results 
of the previous sections, confirming that the effective 
low-energy behavior of a ferrimagnetic chain corresponds 
to a gapless ferromagnet. The approximate ground 
state energy per cell which is obtained from Eq. ( [so]) is 
(l/N)E = -J(l + (<5/9)). As a result of this effective 
description, one can also obtain the ferrimagnetic gapless 
mode dispersion relation (u> ~ k 2 ) by using the disper- 
sion relation of the 5=1/2 Heisenberg ferromagnet in 
Eq. ©. 

In addition to considerations that have to do with the 
ferrimagnetic ground state, we can also gain some in- 
sight about the low-energy excitations of the quantum 
ferrimagnets by constructing another embedding oper- 
ator (Q) which gives the effective Hamiltonian for the 
excited spectrum. The first excited state is obtained by 
considering one dimer in the 5 = 3/2 multiplet Hilbert 
space. In this case the embedding operator of the excited 
dimer is constructed by the base kets of the 5 = 3/2 mul- 
tiplet and the other ones remain as before. The resulting 
Hamiltonian obtained from the QRG process in this case 
is given by 
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ff£)=(-JVJ+|j) + |w(S' 



j'-i ' s "i 



(57) 



where S' = 1/2 and 5" = 3/2. The first term in f57|) con- 
tributes to the gap and the second one is responsible for 
antiferromagnetic interactions. (Therefore the gapped 
excitations exhibit antiferromagnetic interactions). If we 
continue this procedure by considering more dimers in 
the 5 = 3/2 multiplct Hilbert space, other excited states 
are obtained and finally for the highest energy spectrum 
we find 



(58) 



Therefore for sufficiently high energies the physics of the 
(S\ = 1/2, S2 = 1) ferrimagnetic Heisenberg chain corre- 
sponds to a spin 3/2 Heisenberg antiferromagnet which is 
separated by a finite gap from the gapless ferromagnetic 
branch. 

To further substantiate our results we have also cal- 
cualted the (average) magnetization per spin in the unit 
cells and the spin-spin correlation functions. In the 
QRG formalism the ground state |0) is replaced by Q\0'), 
where |0') is the ground state of the effective Hamiltonian 
Hilbert space. Therefore the magnetization, mj, per site 
for each spin is 



to obtain the real-space representation of the effective 
low-energy Hamiltonian explicitly to study the physical 
properties of ferrimagnetic Heisenberg chains. 

The results of our calculation of the correlation func- 
tion between the z-components of the S = 1/2 spins and 
an arbitrary spin on the chain is plotted in Fig. ^. We 
have also plotted the QMC results (S. Brehmer, et.al 
||) for comparison. The QRG results are in complete 
qualitative agreement with the QMC data, but there 
are slight quantitative differences between the two ap- 
proaches. In particular, in the long-wavelength limit the 
QRG approach agrees more closely with the spin-wave 
theory than QMC calculations. However, both QRG and 
QMC yield a correlation lenght which is smaller than a 
lattice spacing. 

In addition to supporting the ideas that were discussed 
in the previous sections, the results in Table. 1 and Fig. |] 
confirm that the QRG method leads to good qualitative 
and quantitative results while maintaining an overall an- 
lytical simplicity. In fact, besides the (1,1/2) that we 
just discussed in detail, we have also studied two other 
cases (Si = 1/2, S 2 = 3/2) and (Si = 1,S 2 = 2) with 
comparable ease. We find that the same general picture 
holds for these two cases as well. Then the effective low- 
energy Hamiltonian for the (Si = 1/2, S2 = 3/2) is given 
by 

H^(l/2, 3/2) = -&NJ - £ Sf > • Sffi, 



mi 



1 N 



(59) 



with I — 1,2 for Si and S2 respectively. (Note that the 
spin- wave sum rule, Mo = mi + m 2 = |Si — S2I, is sat- 
isfied by Eq. (p9|).) We have summerized the results of 
calculations of magnetization, energy gap(A), and the 
ground state energy which has been obtained by using 
the QRG techniques (this section), spin- wave analysis 
(Sec. II and III, f|), DMRG (S. K. Pati et al. g), and 
quantum Monte-Carlo (QMC) calculations (S. Brehmer, 
et al. ||) in Table. 1 for a quantitative comparison. As 
can be seen from Table. 1, the QRG results are closer to 
those found from the spin- wave approach which is, how- 
ever, not as accurate as the DMRG and QMC results. 
It is known p3 that in the QRG method the number 
of kept states and the lenght of blocks is responsible for 
the accuracy of the results due to the boundary effects in 
isolated blocks. However, this loss of accuracy, as com- 
pared to DMRG and QMC, e.g., is compensated for by 
straightforward analytically tractable results which ex- 
plicitly exhibit renormalization group flows towards the 
fixed points. In the example discussed here we were able 



where J, 



(1/2,3/2) 



(60) 

(5/16) JS > and S' 1 ' is a spin 1. For 



the (Si = 1, S2 = 2) case we find 



N 



H%(l,2) = -3NJ-J™Y, S ? Ks & 



(61) 



i=l 



where = (3/4) J 8 > 0. One can use © and @ 

to obtain the ferromagnetic gapless dispersion relations 
(u> r~j k 2 ) for these ferrimagnetics. For example, the ap- 
proximate ground state energy and the magnetization per 
spin are 



and 



±E (l/2, 3/2) = -ij(l + £), 

1 3 
m(-) = -0.25 ; m(-) = 1.25, 
y 2 2 



l£ (l,2) = -3J(l + £), 
m (l) = -0.5 ; m(2) = 1.5. 



(62) 



(63) 



The general conclusion that emerges from the QRG 
calculations that we have performed is that the effec- 
tive low-energy Hamiltonian of a (Si,S 2 ) ferrimagnetic 
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Heisenberg chain is the same as that of a spin S = 
\S\ — S2I ferromagnetic chain with a gapless mode. Our 
QRG study also exhibits the cross-over to antiferromag- 
netic behavior at a finite energy gap which was also ob- 
tained in the previous sections by using semiclassical ap- 
proximations. A more refined analysis of quantum fer- 
rimagnetic spin system within the QRG framework and 
at finite temperature can, in principle, lead to quantita- 
tive predictions in the region where the fcrromagnctic- 
antifcrromagnetic cross-over takes place. In view of the 
complexity of the problem through using the continuous 
renormalization group approach — even when combined 
with semiclassical approximations such as the large-S 1 
limit, e.g. — it would be very instructive to examine these 
questions in future investigations by using QRG tech- 
niques. 



magnetic Heisenberg chain is the same as that of a spin 
S = \Si — S2 1 ferromagnetic chain. The methods and 
applications discussed in this paper have enabled us to 
obtain good quantitative and qualitative results and can 
be extended to study more general quantum ferrimag- 
nets. 
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VII. CONCLUSION 



In this paper we have presented an extensive study 
of quantum ferrimagnets in one, two, and three dimen- 
sions at both zero and finite temperature by using a va- 
riety of techniques. These include the spin-wave theory 
semiclassical approximation, the continuum non-linear 
cr-model type effective field theory, a formulation based 
on the quantum rotor models, and the real space quan- 
tum renormalization group approach. We find a rich 
phase diagram for quantum ferrimagnets with nearest- 
neighbor interactions at T = and T > 0. The effec- 
tive description of these systems considerably simplifies 
in the long- and short- wavelength limits. In the long- 
wavelength limit quantum ferrimagnets behave as ferro- 
magnets with a (quantized) magnetization 15*1 — S2I per 
unit cell, whereas in the short-wavelength limit they ex- 
hibit strong antiferromagnetic features. We also demon- 
strate the existence of an intermediate continuous phase 
transition to a partially polarized magnetized state in 
quantum ferrimagnets. Similarly to the usaul ferro- 
magnets and antiferromagnets, one-dimensional quan- 
tum ferrimagnets do not exhibit long range order at 
non-zero temperatures; however, two-dimensional ferri- 
magnets show a quasi-long range order at non-zero tem- 
perature and a Berezinskii-Kosterlitz-Thouless type of 
phase transition can be seen. Three-dimensional quan- 
tum ferrimagnets behave differently from their one- and 
two-dimensional counterparts and exhibit both ferromag- 
netic and antiferromagnetic long range order, the relative 
strength of which depends on the temperature and the 
ratio S1/S2 of neighboring spins. At T — quantum 
ferrimagnets show long range order at any dimension. 
When the directions of Si and 5 2 are opposite, both 
the magnetization M = \Si — S2I and the staggered 
magnetization M Neel are finite and can be considered 
as order parameters. And finally, We have shown that 
the effective low-energy Hamiltonian of a (Si,^) ferri- 
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FIG. 1. The finite temperature phase diagram for the 
3-dimensional quantum ferrimagncts. 



FIG. 2. The magnetization per unit cell (Mo) as a function 
of the external magnetic field. When h c i < h < h C 2 the 
ferrimagnet is partially polarized. At h — h C 2 the external 
field leads to a saturated polarized (ferromagnetic) state. 



FIG. 3. Block decomposition of the ferrimagnetic chain 
into unit cells. 



FIG. 4. The (Si z S2 Z ) correlation function vs. distance as 
obtained from QRG and QMC. Here Si z = 1/2 and S2z is the 
z-component of spin on an arbitrary site. 



TABLE I. A comparison of the magnetization per site 
(mi), the energy gap (A), and the ground state energy per 
unit cell (Eo/JN) as obtained from QRG, spin- wave (SW) 
theory, DMRG, and QMC calculations. 





QRG 


SW 


DMRG 


QMC 


m(S= ±) 


-0.166 


-0.195 


-0.292 


-0.290 


m(S = l) 


0.666 


0.695 


0.792 


0.790 


A/J 


1 


1 


1.279 


1.767 


Eo/JN 


-1.111 


-1.436 


-1.454 


-1.455 
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